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ABSTRACT 

We present results from our Very Large Telescope large program to study the dynamical 
evolution of local Ultraluminous Infrared Galaxies (ULIRGs) and QSOs. This paper is the second 
in a series presenting the stellar kinematics of 54 ULIRGs, derived from high resolution, long-slit 
H- and if -band spectroscopy. The data presented here, including observations of 17 new targets, 
are mainly focused on sources that have coalesced into a single nucleus. The stellar kinematics, 
extracted from the CO ro-vibrational bandheads in our spectra, indicate that ULIRG remnants 
are dynamically heated systems with a mean dispersion of 161 kms - . The combination of 
kinematic, structural, and photometric properties of the remnants indicate that they mostly 
originate from major mergers and that they result in the formation of systems supported by 
random motions, therefore, elliptical galaxies. The peak of the velocity dispersion distribution 
and the locus of ULIRGs on the fundamental plane of early-type galaxies indicate that the end 
products of ultraluminous mergers are typically moderate-mass ellipticals (of stellar mass ~10 10 - 
10 11 Mq). Converting the host dispersion into black hole mass with the aid of the Mbh — cr 
relation yields black hole mass estimates of the order 10 7 - 10 s Mq and high accretion rates with 
Eddington efficiencies often > 0.5. 



Subject headings: galaxies: formation 
ISM: kinematics and dynamics — 

Introduction 



galaxies: kinematics and dynamics — infrared: galaxies 



Galaxy mergers, which have a frequency in- 
creasing with redshift (e.g., Toomre 1977; Kauff- 
mann & White 1993; Le Fevre et al. 2000), are 
considered a key mechanism in driving galaxy evo- 
lution. In the local Universe, the ultraluminous 
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infrared galaxies (ULIRGs) are excellent laborato- 
ries for studying violent merging events. ULIRGs 
are mergers of gas-rich galaxies observed during 
strong starburst events; these episodes have dura- 
tions <10 8 yrs (e.g. Canalizo k Stockton 2001; 
Mihos k Hernquist 1996) and are brief compared 
to the baryonic matter merging process (~10 9 yrs; 
e.g. Barnes 1992; Hernquist 1993). The star- 
burst emission, often combined with emission orig- 
inating from an active galactic nucleus (AGN), 
gives rise to infrared (IR) luminosities greater than 
lO 12 ^© (e.g., Sanders & Mirabel 1996; Lonsdale 
et al. 2006 and references therein); these lumi- 
nosities are comparable to the bolometric lumi- 
nosities of QSOs. 

A plethora of studies indicates that ULIRGs 
transform (gas-rich) spiral (S) galaxies into ellipti- 



1 



cals (Es) through merger-induced dissipative col- 
lapse (Kormendy & Sanders 1992; Mihos & Hern- 
quist 1996; Barnes & Hernquist 1996). The large 
molecular gas concentrations in the central kpc re- 
gions of ULIRGs (e.g. Downes & Solomon 1998; 
Bryant & Scoville 1999) have densities compara- 
ble to stellar densities in ellipticals. The light pro- 
files of ULIRG remnants often follow a r 1 / 4 law 
(e.g., Wright et al. 1990; Stanford & Bushouse 
1991; James et al. 1999; Scoville et al. 2000; 
Rothberg & Joseph 2004). Kim et al. (2002) 
and Vcilleux et al. (2002) established the so- 
lidity of this result by proving that the majority 
(73%) of all single-nucleus ULIRGs in a sample 
of 118 sources with 60 /um flux greater than 1 
Jy (hereafter the 1 Jy sample; Kim & Sanders 
1998) have elliptical-like light profiles. Genzel 
et al. (2001) and Tacconi et al. (2002) have 
made high-resolution near-infrared (NIR) spectro- 
scopic measurements of the stellar dynamics of 
small samples consisting mostly of fully-merged 
ULIRGs. They conclude that ULIRGs resem- 
ble intermediate mass ellipticals/lenticulars with 
moderate rotation, in terms of their velocity dis- 
persion distribution, their location in the funda- 
mental plane (FP; e.g., Djorgovski & Davis 1987; 
Dressier et al. 1987) and their distribution of the 
rotation/ velocity dispersion ratio. Together, these 
results suggest that ULIRGs form moderate mass 
ellipticals (of stellar mass ^lO 11 M Q ). 

One way to investigate the physical details and 
the evolution of ULIRGs is to determine the kine- 
matic and structural properties of the merging 
galaxies in different merger phases. We there- 
fore conducted a European Southern Observatory 
(ESO) large program 2 that traces the host dynam- 
ics of a large sample of ULIRGs spanning wide 
ranges of merger phase and infrared luminosity 
through NIR spectroscopy. Our study extends the 
previous work of Genzel et al. (2001) and Tac- 
coni et al. (2002). The enlarged sample comprises 
ULIRGs at wider ranges of merger phase and lu- 
minosity, and higher redshifts than those of Gen- 
zel et al. (2001) and Tacconi et al. (2002). The 
new observations of 38 sources increase the origi- 
nal sample by more than a factor of 3, enhancing 
our statistics and enabling us to study the prop- 
erties of galaxy mergers as a function of time. In 
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Dasyra et al. (2006, hereafter Paper I) we ana- 
lyzed those ULIRGs that are in a merger phase 
later than the first encounter but prior to nuclear 
coalescence, and hence, show more than one nu- 
cleus in the NIR acquisition images. In this paper 
we present results mainly from those ULIRGs that 
have coalesced and show a single nucleus in our 
images, the so-called remnants. We compare the 
stellar kinematic properties of binary ULIRGs and 
ULIRG remnants to look for traces of dynamical 
evolution in ultraluminous mergers. 

Some of the ULIRGs presented in this study 
may in fact be binary sources very close to co- 
alescence that have small projected nuclear sep- 
arations; such sources cannot always be resolved 
or kinematically disentangled due to instrumental 
angular resolution constraints. At redshifts typi- 
cal for the sources in our sample, the angular res- 
olution achieved implies that any unresolved sys- 
tems will have nuclear separations smaller than 1.5 
kpc. Merger simulations (e.g. Mihos 1999; Mihos 
2000) have shown that by the time the individual 
nuclei are separated by < 1 kpc, the low moments 
of the stellar kinematics (rotational velocity and 
velocity dispersion) have almost reached their re- 
laxation values. Therefore, the dynamical proper- 
ties of all the sources we classify as remnants are 
representative of those at dynamical equilibrium. 

This paper is arranged as follows. We briefly 
summarize the observations and data reduction 
methods and present the pre- and post- coales- 
cence ULIRG samples in § 2. After studying 
whether the kinematic properties of the merging 
galaxies evolve during the ultraluminous merger 
phases in § 3, we investigate the origins and the 
potential end products of ultraluminous mergers 
in § 4 and § 5 respectively. We then focus on 
the black holes hosted by ULIRGs: an analysis of 
the evolution of the Mbh- o relation during the 
merger is followed by a discussion of the nuclear 
activity implied by our data in § 6. Finally a sum- 
mary is presented in § 7. 

2. Observations And Data Reduction 
2.1. Data Acquisition And Analysis 

We combine the sources in our Very Large Tele- 
scope (VLT) program with the sources presented 
in Genzel et al. (2001) and Tacconi et al. (2002) 
to compile a sample that comprises 54 ULIRGs 
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and 12 QSOs. The stellar kinematics of the QSOs 
will be presented in a forthcoming paper (Dasyra 
et al. 2006, in preparation). In this paper, we 
perform a study of the stellar kinematics in post- 
coalescence ultraluminous mergers (ULIRG rem- 
nants), after presenting new data for 17 ULIRGs 
and summarizing pre-existing data for 13 mainly 
post-coalescence sources. Of the 54 ULIRGs ob- 
served in total, 30 are confirmed to have two nu- 
clei 3 , and 1 (IRAS 00199-7426: Due et al. 1997; 
Paper I) may be a multiple-interaction system. A 
detailed description of the criteria set to select 
these 54 sources from the 1 Jy catalog is given 
in Paper I. 

Our new high-resolution, long-slit spectroscopic 
data were obtained using the ISAAC spectrome- 
ter (Moorwood et al. 1998) mounted on the Antu 
telescope unit of the VLT. The observations were 
made in the H and K bands, with an instrument 
resolution of A/<5A=5100 and A/<5A=4400 respec- 
tively. The slit length was 120" and the slit width 
0'.'6; the detector scale was 0'.'146 per pixel. The 
total exposure times and the slit position angles 
for each source are presented in Table 1. If the 
position angle of the major axis of rotation could 
not be identified for the remnants (e.g. from the 
elongation of the stellar disk), the slits were typ- 
ically placed at 0° and 90°. The redshift range 
of all sources in our sample is 0.018 < z < 0.268 
(see Table 1). No sources are observed for the red- 
shift sub-range 0.163 < z < 0.199, since the CO 
bandheads are then shifted into wavebands of high 
atmospheric absorption. 

As in Paper I, we derive the structural parame- 
ters of the sources that do not have high-resolution 
NIR imaging (e.g. Hubble Space Telescope NIC- 
MOS) data by fitting ellipses to our 10-second- 
long acquisition images. The acquisition images 
have also been obtained with the ISAAC detector, 
in its H— band imaging mode, which corresponds 
to a scale of 0'.'147 per pixel. We fit the ellipses 
with the aid of the SExtractor package (Bertin & 
Arnouts 1996) and we present the half-light radius 
R e g 4 , the ellipticity, and the angle <\> a between 
the major axis of rotation and the position angle 
of the first slit in Table 2 for each source. To con- 
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vert all angular distances into linear sizes we use 
a Ho=70 km s _1 Mpc -1 , fl m =0.3, f2totai=l cos- 
mology. For some of the sources in our sample, 
a measurement of the if-band effective radius is 
available in the literature (see Table 2). In these 
cases, we use the average NIR i? e ff value in all com- 
putations and diagrams that follow. The imaging 
analyses of Surace & Sanders (1999) and Scovillc 
et al. (2000) have shown that the difference be- 
tween the measurement of R c s in ULIRGs in the 
H and in the K band is small and non-systematic. 
Therefore, averaging these NIR i? c ff values yields 
a measurement of higher accuracy. 

The extraction of the stellar central velocity 
dispersion a and rotational velocity V TOt from the 
spectra follows the method presented in Paper I. 
It is performed using the Fourier quotient tech- 
nique (Bender 1990); this method provides the 
intrinsic line-of-sight (LOS) velocity profile along 
any given aperture. To this we fit a high-order 
Gaussian (linear combination of a Gaussian and 
a second order polynomial) to determine the av- 
erage LOS radial velocity and velocity dispersion. 
The fit is performed to each bandhead individu- 
ally and the errors are equal to the standard devi- 
ation of all measurements performed. The central 
aperture spectrum of each source, combined over 
the two slits and shifted to restframe, is displayed 
in Fig. 1. The stellar template that fits best the 
post-coalescence ULIRG spectra is an M0III gi- 
ant (HD 25472; presented in Gcnzcl et al. 2001). 
In general, template mismatch affects the accu- 
racy to which the stellar velocities are measured; 
we found that the difference in the velocity disper- 
sion between the (best-fitting) M0III giant and the 
Mil supergiant (HD 99817; Genzel et al. 2001) is 
~15 kms -1 . The M0III template is overplotted 
with a dashed line in Fig. 1, after being convolved 
with the Gaussian that best simulates the LOS 
broadening function. 

In all of this paper's analyses we deem aper- 
ture effects to be negligible for the measurement 
of the stellar velocity dispersion. According to the 
merger simulations of Bendo & Barnes (2000), 
the velocity dispersion in (relaxed) remnants may 
vary at most by 10% for the apertures considered 
here (up to ~2i? e ff) which are selected to maxi- 
mize the signal-to-noise ratio. Resolution effects 
further reduce differences in the velocity disper- 
sion measurement between various apertures. 
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The LOS rotational velocity, U ro t(LOS), is mea- 
sured from apertures that exclude the center of 
the galaxy. The center and the annular width of 
the apertures used are tabulated in Table 3. For 
each slit, the observed rotational velocity V rot {obs) 
can be derived from the LOS rotational velocity 
by correcting for the angular deviation of the slit 
from the major axis of rotation as 

V Iot (obs) = V Iot [LOS}/cos[<f> a }. (1) 

When rotation is observed along both slits, i.e. 
when none of the slits is very close to being per- 
pendicular to the major axis of rotation, we aver- 
age the results of the two slits. The value of the 
observed rotational velocity is given in Table 3 and 
is related to the actual rotational velocity V TOt by 

l/ rot (obs) = V Iot sin(i). (2) 

In this paper we do not compute the inclination 
i from the ellipticity of each source (as we did in 
Paper I); this conversion is very uncertain for the 
ULIRG remnants since their stellar disks are dy- 
namically hot (due to the advanced phase of the 
merger). Instead, we statistically determine the 
mean inclination that needs to be applied to the 
sample. For this purpose, we use the weighted 
mean value of sin(i) observed for disk-like galax- 
ies on the sky. Since the probability p[i] of finding 
a galaxy at an inclination i (in the range [0°,90°]) 
scales with sin(i) (Collin et al. 2006), the mean 
value of p[i]sin(i) is 7r/4. The stellar kinematic 
results (central velocity dispersion and rotational 
velocity) can be found in Table 3. 

2.2. Pre- and Post-coalescence Sample 
Classification 

Prior to performing statistics on the kine- 
matic parameters of the pre- and post- coalescence 
ULIRG samples, we descibe the criterion we use 
to compile the two samples; it consists of a nuclear 
separation threshold beyond which binary sources 
can be considered relaxed. According to simula- 
tions, this nuclear separation should be roughly 1 
kpc (Mihos 2000). The actual value of the thresh- 
old we set depends on the resolution of the images 
we use to search for the presence of secondary 
components. 

We compute the resolution of our images, which 
depends on the detector scale and the seeing, by 



averaging the FWHM of bright stars in each field. 
Our mean resolution corresponds to a FWHM of 
4.3 pixels and enables us to resolve individual 
sources separated by 1 kpc up to z=0.080. If we 
set the nuclear separation threshold to 1.5 kpc, we 
can distinguish individual sources up to 2=0.118. 
The corresponding numbers of sources that are 
verified to satisfy the coalescence criterion are 10 
and 23, respectively. We adopt a nuclear separa- 
tion cutoff of 1.5 kpc. For this nuclear separation, 
the low velocity moments are still close to their re- 
laxation values according to the models of Mihos 
(2000) and Naab et al. (2006b). 

Of the 7 sources at z > 0.118, 5 have 
been observed with the NICMOS camera on- 
board the Hubble Space Telescope by Veilleux 
et al. (2006) and Scoville et al. (2000). Ac- 
cording to these authors, none of IRAS 00397- 
1312, IRAS 01572+0009, IRAS 04313-1649, IRAS 
09039+0503, and IRAS 14070+0525 has a sec- 
ondary object at distances <1.5 kpc that can be 
unambiguously characterized nucleus. For 

one of the remaining 2 sources, IRAS 23578-5307, 
the large elongation and the tidal tails indicate the 
possible presence of two components (see Fig. 2; 
left panel). To investigate this possibility, we 
deconvolve our acquisition image with its point- 
spread- function (PSF) which has a FWHM of 4.1 
pixels. For this task, we use the LUCY algorithm 
of IRAF (Lucy 1974; Richardson 1972). The de- 
convolution increases the resolution of the image 
by a factor of ^2. The resulting image (see Fig. 2; 
right panel) confirms the presence of a second nu- 
cleus and constrains its projected separation to 
be 1.4 ± 0.4 kpc. Unlike that of IRAS 23578- 
5307, the redshift of IRAS 11223-1244 is too high 
to enable us to derive any conclusions using im- 
age deconvolution techniques. IRAS 11223-1244 is 
therefore excluded from all statistics that follow. 
The same applies to the possibly multiple merger 
IRAS 00199-7426. 

A summary of the classification of all sources 
in our sample can be found in Table 1. We 
have a total sample of 29 ULIRG remnants, 
which consists of 21 single- nucleus or unresolved 
sources, and 8 confirmed close binaries that 
have nuclear separations less than or equal to 
the selected threshold. These sources namely 
are IRAS 00091-0738, IRAS 09111-1007, IRAS 
14378-3651, IRAS 15250+3609, IRAS 23578-5307, 
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Mrk 273, Arp 220 and NGC 6240 (see the Ap- 
pendix). The pre-coalescence binary sample com- 
prises 23 pairs of sources; of those, 40 individ- 
ual components have velocity dispersion measure- 
ments. 

3. Traces of Evolution In The Stellar Kine- 
matics 

Taking into account the classification scheme 
of § 2.2, we find that the a distribution of the 
29 ULIRG remnants has a mean of 161 kms -1 
and a median of 150 kms -1 . The uncertainty in 
this mean, which is known as its standard error, 
is 8 kms -1 since the standard deviation equals 
42 kms -1 . For the 40 pre-coalescence ULIRG nu- 
clei, we find that the mean is 142 kms -1 , with a 
standard error of 3 and a standard deviation of 21 
kms -1 , and the median is 145 kms -1 . The pre- 
and post- coalescence dispersion distributions are 
shown in Fig. 3 (left panel). While the difference 
in the mean of the two distributions is small and 
the modes are the same, the remnant velocity dis- 
persion distribution has a larger variance and a 
tail at the high-a end, the statistical significance 
of which needs to be quantified. 

To investigate whether the two distributions are 
independent, we begin by using the Kolmogorov- 
Smirnov (KS) test, which makes no assumptions 
about the shape of the distributions under exam- 
ination. We find that the maximum deviation D 
between the cumulative distribution function of 
the two ULIRG populations is 0.285, which for the 
number of sources we observed, corresponds to a 
probability of 89.4% that the populations are inde- 
pendent. This probability is below the widely ac- 
cepted significance levels. However, a well-known 
disadvantage of the KS test is that it is sensitive 
to small number statistics. To address whether 
a better handling of small number statistics may 
alter our conclusions, we created Monte Carlo sim- 
ulations. 

The Monte Carlo code begins with the assump- 
tion that the two distributions are drawn from the 
same parent population. We simulated the par- 
ent distribution by spline interpolating between 
the observed a distribution of the combined pre- 
and post- coalescence ULIRGs (see Fig. 3, right 
panel). From the parent distribution, which con- 
sisted of 10 6 points, we selected random points 



to generate two artificial progenitor/remnant sub- 
samples. Each generated subsample had a num- 
ber of elements equal to that of the corresponding 
real population. We repeated this procedure for 
10000 iterations. In each iteration, we computed 
the difference in the mean of the two generated 
samples and the difference in the variance be- 
tween each real population and the respective 
generated sample. At the end of all iterations, 
we counted how many times the difference in the 
mean of the two generated samples was equal to 
or greater than that measured, 19 kms -1 . We 
found that it corresponded to a probability of 
1.2%. Similarly, the probability p that the vari- 
ances of both generated samples (c| e n,i an d a gen,2 
respectively) are greater than those measured 
(or p[cr g on4 > 21kms~ 1 ] *p[er g on,2 > 42kms -1 ]) 
is 5.0%. In other words, the different means imply 
a 98.8% probability that the two distributions are 
independent, while the different variances yield a 
95.0% probability of independence. Since these 
probabilities correspond to ~2.5 and ~2.0 sigma 
respectively, the difference in the distributions is 
at or above the widely used significance levels. 
This result favors the hypothesis that the appar- 
ent kinematic evolution is real rather than an 
artifact of limited-number statistics. 

Despite the large range of confidence levels in- 
dicated by the statistical tests, their conclusion 
can be summarized as follows: it is possible, al- 
though uncertain, that the observed increase in 
sigma reflects the dynamical heating of the merg- 
ing galaxies. The importance of this result lies in 
that observations are able to constrain the dynam- 
ical heating that models predict for the ultralumi- 
nous phases of gas-rich mergers. Theory (Mihos 
& Hernquist 1996; Mihos 1999; Springel et al. 
2005) and observations (Mihos & Bothun 1998; 
Murphy et al. 2001; Veilleux et al. 2002) show 
that strong starbursts (and associated ULIRG 
phases) in gas-rich mergers typically occur be- 
tween first encounter and nuclear coalescence, al- 
tough in individual cases, their actual occurence 
may significantly vary in terms of time. On a sta- 
tistical basis, a binary ULIRG sample is expected 
to be roughly mid-way from first encounter to co- 
alescence; a ULIRG remnant sample is probably 
close to coalescence. For these merger phases, the 
theoretically predicted increase of the stellar ve- 
locity dispersion is small (Naab et al. 2006b); 
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therefore observationally tracing it is very impor- 
tant. However, prior to comparing observations 
with models, further data are needed to quantify 
at 3 Gaussian-sigma levels the difference in a be- 
fore and after nuclear coalescence. 

Another possibility that needs to be investi- 
gated prior to comparing the evolution in the stel- 
lar kinematics between observations and models is 
whether the double- and single- nucleus ULIRGs 
may originate from mergers of galaxies of some- 
what different mass contents: the ULIRGs that 
are prior to nuclear coalescence may be more gas- 
rich (and possibly more massive) than the rem- 
nants, allowing them to have an ultraluminous 
phase of comparable IR output at earlier merger 
phases (Mihos & Hcrnquist 1996). In this case, 
the observed dynamical heating would be less than 
what would be measured if we were able to ob- 
serve a specific galaxy pair from the beginning 
until the end of the merging process. To inves- 
tigate whether the starburst activity is triggered 
under similar conditions for pre- and post- coales- 
cence ULIRGs, we need to quantify and compare 
the gas-mass content of both samples. Molecu- 
lar gas mass measurements obtained by Gao & 
Solomon (1999) in local Luminous Infrared Galax- 
ies (LIRGs; sources of 10 n L© < L m < 1O 12 L ) 
and ULIRGs indicate that a correlation between 
nuclear separation and gas fraction is observed in 
LIRGs but not in ULIRGs. To properly address 
the question, molccular-gas-content observations 
need to be performed on ULIRG samples with 
sizes comparable to those of the dynamical stud- 
ies. 

The mean value of the rotational velocity is 62 
(and the standard deviation is 59) kms -1 for the 
ULIRG remnants, increasing to 79 kms -1 when 
inclination effects are statistically accounted for. 
The observed stellar rotational velocity of each 
individual source and its ratio over the velocity 
dispersion is presented in Table 3. The mean 
V TO t (obs)/<7 ratio of the remnants is 0.36, increas- 
ing to V rot /a—0A6 when we apply the statisti- 
cal inclination correction to the rotational veloc- 
ity. The individual progenitors of the binary sys- 
tems have a mean rotational velocity of 59 (with 
a standard deviation of 38) kms -1 or 105 (with 
a standard deviation of 96) kms -1 when the in- 
clination effects are corrected from the ellipticity 
of each progenitor. The corresponding observed 



and inclination-corrected V ro t /a ratios of the pro- 
genitors are 0.42 and 0.76. The somewhat lower 
V mt I & ratio of the remnants can only be attributed 
to their higher (than the progenitors) value of u, 
since the difference between the pre- and post- co- 
alescence rotational velocity is insignificant. In 
these calculations we have not attempted to cor- 
rect the central velocity dispersion for inclination 
effects. 

4. Origin Of The ULIRG Remnants 

We infer the initial conditions of the mergers 
that lead to ultraluminous IR activity by compar- 
ing the kinematic properties of ULIRG remnants 
with those predicted by simulations in the litera- 
ture. According to various authors (e.g. Bendo & 
Barnes 2000; Naab & Burkert 2003) the remnant 
V4ot(obs)/cr ratio depends upon the initial mass 
ratio of the merging galaxies. In the gas-free, N- 
body simulations of binary mergers performed by 
Naab & Burkert (2003) for several mass ratios and 
orientations, the major mergers produced slowly 
rotating remnants. Those authors suggested that 
the U ro t(obs)/(J ratio is ~ 0.2 for 1:1 and ~ 0.4 
for 2:1 merger remnants, while it reaches higher 
values (0.8) when the remnants originate from mi- 
nor (4:1) mergers. Our remnants agree best with 
a 1:1 and 2:1 merger origin, also agreeing with the 
directly measured progenitor mass ratios of Pa- 
per I and confirming that ULIRGs are representa- 
tive of the most violent local mergers. The results 
are similar when comparing to those of Naab et 
al. (2006b), who have presented similar contours 
for gas-rich, dissipative simulations of 1:1 and 3:1 
mass ratios. 

Naab & Burkert (2003), Burkert & Naab 
(2005), and Naab et al. (2006b) have also shown 
a connection between the V rot (obs)/cr ratio and 
the remnant ellipticity at the effective radius for 
several progenitor mass ratios. We overplot our 
results with those of the merger simulations in 
Fig. 4, where our ULIRGs are given in triangles. 
Each panel corresponds to mergers of different 
progenitor mass ratios (denoted at the upper-left 
corner of each panel). In bold, solid, and dotted 
contours we show the 90%, 70%, and 50% proba- 
bility of finding a merger remnant of each category 
in the enclosed region. These contours correspond 
to models that do not include gas (Naab & Burk- 
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ert 2003). For the dissipative merger remnants of 
Naab et al. (2006b), the 90 % and 50% probabili- 
ties are shown as dark- and light-gray shaded areas 
respectively. A comparison of the predicted with 
the measured remnant dynamics indicates that 
ULIRGs are more likely to be formed by major 
mergers (see Paper I). We also find that ULIRGs 
have a narrower mass-ratio origin than the merger 
remnants in Rothberg & Joseph (2006b), who 
have carried out a similar analysis on a sample 
of sources that are selected visually to have per- 
turbed morphologies (see also § 5.3). 

The addition of gas to the simulations of Naab 
& Burkcrt (2003) did not significantly influence 
the position of remnants on the V TOt (obs)/a vs. 
e diagram (Burkert & Naab 2005; Naab et al. 
2006b). However, the addition of star formation 
may result in an increase of the theoretically pre- 
dicted Kot(obs)/<7 ratio and ellipticity, since new 
stars will be created upon a dynamically-cold disk, 
formed by gas that is not consumed during the 
merger (Barnes 2002; Springel 2000). This pos- 
sible increase would shift the probability contours 
to higher values on both axes further supporting 
our findings that ULIRGs are produced mainly by 
mergers of galaxies of comparable masses. On the 
other hand, the position of ULIRGs on this dia- 
gram could be shifted upwards along the V TOt /a 
axis due to aperture correction effects. The V mt 
values used are not necessarily measured from 
apertures where the rotation curve has reached its 
flat part (van Albada et al. 1985). However, 
the facts that V rot has been measured from aper- 
tures between 0.5 and 2 R e s and that the elliptic- 
ity at the effective radius does not depend upon 
the apertures used to measure V TOt are good indi- 
cations that such corrections will probably not be 
sufficiently large to alter our conclusions. 

5. End Products Of ULIRGs 

5.1. The Masses Of Ultraluminous Merger 
Remnants 

A direct way to investigate the nature of the 
end products of ultraluminous mergers is to cal- 
culate the mass m of their remnants. According 
to Bender et al. (1992) the bulge mass is related 
to the stellar kinematics as 

m = c 2 af 00 R cS , (3) 



where cioo is the projected central velocity disper- 
sion in units of 100 km s _1 , i? e ff is in kpc, and to is 
in 10 10 Mq. The geometrical factor C2 depends on 
the distribution of matter in the galaxy. Following 
Tacconi et al. (2002), we adopt C2 = 1.4, as ap- 
propriate for a constant ra/L King model whose 
tidal-to-core radius ratio is 50, midway between 
those of dwarf and giant ellipticals. By combining 
the above, the total dynamical mass is computed 
from 

to = 4.7 x 10 5 (3(J 2 + V mt 2 )R cS , (4) 

where a and V Iot are now in units of kms -1 i? e ff 
is in kpc, and to is in Mq. 

We use the mean values of the stellar disper- 
sion and inclination-corrected rotational velocity 
from § 3. For the objects that do not have a 
V TOt measurement (see Table 3), we use their dis- 
persion in combination with the mean inclination- 
corrected Vrot/c ratio (0.46) to infer it. We find 
that the mean dynamical mass of the remnants is 
8.91 x 10 10 M Q (in good agreement with Tacconi et 
al. 2002), suggesting that ultraluminous activity 
mainly originates from mergers of sub-m* galaxies, 
for to* = 1.4 x 10 11 Mq (Genzel et al. 2001 and 
references therein). For different values of to* (e.g 
that of Bell et al. 2003 adapted to our cosmol- 
ogy, 8.3 x 10 10 M Q ) the sub-TO* characterization 
of ULIRGs may change to <~ to* . Still, in the local 
Universe, the disk progenitors of ULIRGs do not 
need to be as massive as e.g. the Milky Way. 

ULIRG masses derived from Ha emission-line 
dispersions are also sub-m* (Colina et al. 2005). 
Tracing the ionized gas by the [Fell] emission lines 
that appear in our spectra, we find that while the 
mean dispersion of the gas is similar to that of 
the stars, in individual cases gas and stellar dis- 
persions may significantly deviate due to gas out- 
flows. Gas kinematics in these ULIRGs will be 
presented in a forthcoming paper (Tacconi et al. 
2006, in preparation). 

The result that ULIRGs are sub-m* to to* 
galaxies does not contradict findings from imag- 
ing studies that ULIRGs have NIR luminosities 
greater than L* (e.g. Sanders et al. 2000; Col- 
ina et al. 2001; Veilleux et al. 2006). ULIRGs 
can simultaneously be sub-TO* and (1-2) xL* (in 
the .ff-band) mergers since they are selected at 
the peak of their starburst, at which time their 
luminosity-to-mass ratio rises (see also Tacconi et 
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al. 2002; Rupke et al. 2002; 2005a; 2005b). 

5.2. ULIRGs And The Fundamental Plane 
Of Early-type Galaxies 

The V rot I cr ratios of the ULIRG remnants indi- 
cate that mergers of ultraluminous infrared output 
lead to random-motion-dominated systems with 
a non-negligible rotational component. To inves- 
tigate what type of ellipticals ultraluminous IR 
mergers form, Genzel et al. (2001) and Tacconi 
et al. (2002) placed our initial sample of ULIRGs 
on the fundamental plane of early-type galaxies 
(Djorgovski & Davis 1987; Dressier et al. 1987) 
that relates the velocity dispersion, eifective ra- 
dius and surface brightness of these sources. Gen- 
zel et al. (2001) and Tacconi ct al. (2002) con- 
cluded that the remnants resemble moderate mass 
ellipticals (of stellar mass ~10 n M Q ). Our new 
data increase the number of objects in the ULIRG 
luminosity range to be compared with early-type 
galaxies. 

The R c s - a projection of the plane that we 
construct from our data is shown in Fig. 5. Data 
for early-type galaxies are taken from Bender et al. 
(1992), Faber et al. (1997), and Pahre (1999). 
Giant boxy ellipticals (squares) occupy the upper- 
right corner of the FP projection, while disky, 
moderate-mass ellipticals (circles) are located at 
the center. On the left panel of Fig. 5 we overplot 
the ULIRG remnants in triangles; 29 are from this 
study and 2 from Rothberg & Joseph (2006a), 
UGC 5101 5 and AM 2246-490 The location of 
ULIRGs on the plane agrees very well with that 
of moderate-mass ellipticals. On the right panel of 
Fig. 5, we show where remnants of other popula- 
tions are located. Local LIRGs (filled diamonds) 
lie from the lower-left corner of the FP projection 
up to the locus of giant Es (data from Shier & Fis- 
cher 1998, James et al. 1999, Rothberg & Joseph 
2006a, and Hinz & Rieke 2006 6 . The position of 
ULIRGs on this FP projection is more tightly re- 
stricted than that of LIRGs; this result implies 
that ULIRGs have a narrower range of intrinsic 

5 The H- and _R"-band effective radii of this source equal 0.27 
kpc and 1.42 kpc respectively (Scoville et al. 2000; Roth- 
berg & Joseph 2004). 

6 Whenever several velocity dispersion and effective radius 
measurements exist in the literature for sources that are 
not included in our sample, we use the average value of 
those measurements. 



dynamical properties than LIRGs. Other visually- 
selected merger remnants (i.e. sources with per- 
turbed morphology but no IR excess; Rothberg 
& Joseph 2006a) are shown as open-crossed di- 
amonds. That these merger remnants lie closer 
than both the ULIRGs and the LIRGs to the lo- 
cus of giant Es possibly reflects a different merger 
origin (see § 5.3). 

In the 3-dimensional view of the plane, ULIRGs 
are known to deviate from the position of early- 
type galaxies along the surface brightness axis due 
to extinction and population effects (e.g., Genzel 
et al. 2001). Due to the strong starbursts ULIRGs 
undergo, they have a significant population of 
newly formed asymptotic-giant-branch stars, red 
supergiants and giants. Thus, their NIR surface 
brightnesses are higher than those of quiescent el- 
lipticals (Pahre 1999; Veilleux et al. 2002). Prior 
to comparing ULIRGs and Es, the removal of light 
originating from starburst components is therefore 
instructive. Since a significant fraction of the star- 
burst emission is nucleated, we opt to remove the 
central PSF (Veilleux et al. 2006) simultaneously 
removing any AGN-originating emission. 

In the _ff-band, the nuclear PSF removal has 
already been performed by Surace & Sanders 
(1999), Colina et al. (2001) and Veilleux et al. 
(2006) for most of the sources in our sample. For 
the sources taken from Scoville et al. (2000), we 
use the mean if-band ratio between the luminos- 
ity of the PSF-subtracted galaxy and the total 
luminosity, which equals 0.64 (Surace & Sanders 
1999; Colina et al. 2001; Veilleux ct al. 2006). 
The total magnitudes 7 of approximately half of 
the sources we use to construct the 3-dimensional 
view of the iVband fundamental plane arc from 
Kim et al. (2002). For the remaining sources, the 
data are from Due et al. (1997), Surace & Sanders 
(1999), Rigopoulou ct al. (1999), and Scoville et 
al. (2000). In the K band, the value of the ra- 
tio between the luminosity of the PSF-subtracted 
galaxy and the total luminosity is 0.75 times that 
in the H band (Surace & Sanders 1999) i.e., 0.48. 
We compute the mean surface brightness (within 
the effective radius) of each galaxy < yu ga i >, by 
adding 2.blog(2-KR^, s ) mags (where R e g is in ") to 
its PSF-subtracted magnitude in both bands. In 

7 The total magnitude here includes flux integrated up to the 
galaxy's truncation radius. 
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Table 2 wc tabulate the H and the if-band value 
of < /igai > for each source. 

In Fig. 6 (upper panel) we show the 3-dimensional 
view of the fundamental plane in the .ff-band for 
the 21 remnants of Table 2 with a < /J, ga ,i(H — 
band) > measurement. For viewing clarity, we 
also plot in Fig. 6 (lower panel) the fundamental 
plane in one of the representations introduced in 
Pahre (1999). iJ-band data for (mainly cluster) 
Es are from Zibetti et al. (2002) and references 
therein. The ellipticals are plotted as open circles 8 
and the ULIRGs as triangles. In the -ff-band, the 
individual PSF subtraction brings ULIRGs close 
to the fundamental plane, indicating that ULIRG 
remnants will resemble Es once their on-going 
starbursts cease. That ULIRGs are on average 
somewhat brighter than Es is justified by the fact 
that a (diffuse) part of their starburst emission 
remains unremoved after the PSF subtraction. In 
Fig. 7 we show the if-band 3-dimensional view 
and Pahre (1999) visualization of the plane for 25 
remnants in this study with a < a ga \ (K — band) > 
measurement and two sources from Rothberg & 
Joseph (2006a). In the if-band, ULIRGs lie 
again close to moderate-mass Es, after the statis- 
tical PSF removal. Along the photometric axis of 
the plane, PSF-subtracted ULIRGs are also close 
to LIRGs and other merger remnants. 

We have not attempted to correct the photo- 
metric data of any wavelength for extinction ef- 
fects; the value of the usually-assumed equivalent 
screen correction (Scoville et al. 2000; Genzel et 
al. 2001; Tacconi et al. 2002) is uncertain and 
changes after removal of the nuclear PSF. The ap- 
plication of an extinction correction will make the 
surface brightness of ULIRGs brighter and require 
more fading for them to evolve onto the fundamen- 
tal plane. No k-corrections have been applied since 
they are negligible for the majority of the sources 
in our sample. 

5.3. Discussion: The role of ULIRGs in 
the formation of local Es 

Our main conclusion from the fundamental 
plane analysis is that ULIRGs typically lie on 
the locus of local moderate-mass ellipticals whose 



To enhance the clarity of the diagram we have excluded the 
dwarf Es, since they do not fall on the plane of all other 
Es. 



stellar mass is 10 10 — 10 11 M Q . To further inves- 
tigate the type of elliptical galaxy that ULIRGs 
form at the highest rates, we compare the number 
densities of local ULIRG remnants and ellipticals 
as a function of their stellar velocity dispersion. 

In Fig. 8, we plot the source number density 
per stellar dispersion in the range between 50 
and 300 km s" 1 ; ULIRGs are plotted as a his- 
togram and local SDSS ellipticals (0.01 < z < 0.3; 
Bernardi et al. 2003) as a solid line. For the ellip- 
ticals, n(a) is computed from the velocity disper- 
sion function of Sheth et al. (2003). The n(a) dis- 
tribution of ULIRGs is calculated by multiplying 
the % fraction of ULIRG remnants per a bin with 
the volume density of ULIRGs (2.5 x 10~ 7 Mpc~ 3 , 
for our cosmology) from Sanders et al. (2003). 
To facilitate the comparison between the two n(a) 
distributions, wc apply a normalization factor of 
7 x 10 3 to the ULIRG histogram so that its mean 
has the same number density as that of the SDSS 
ellipticals. Physically, this normalization factor is 
related to the ratio between the time over which 
ellipticals have been formed and the lifetime of 
a single ultraluminous burst. Its value is only a 
rough estimate since it depends on the complete- 
ness of the 1 Jy and the SDSS (sub-)samples and 
the fact that local ultraluminous mergers are not 
the only mechanism to produce elliptical galaxies. 

The n(a) distribution of SDSS elliptical galax- 
ies has a mean dispersion similar to that of 
ULIRGs, 157 (with a standard deviation of 54) 
kms -1 . In other words, the descendants of 
ULIRGs will resemble the ellipticals that are most 
common in the local Universe. Local ULIRGs 
seem to form at highest rates sources of character- 
istic <7 between 130 and 160 kms -1 . That local 
Es are mainly produced by mergers of spirals of 
specific luminosities (or masses) is also found in 
the recent simulations of Kaviraz et al. (2006). 

At low dispersions, the ratio between n(a) and 
n(< a >) seems to be lower for ULIRGs than for 
ellipticals. This deviation can probably be at- 
tributed to the fact that mergers of gas-rich galax- 
ies below a certain mass threshold do not possess 
enough gas to undergo a ULIRG phase. 

Deviations between the two distributions also 
exist at the high-dispersion end. For most bins at 
the high-dispersion end, the ULIRG distribution 
seems again underpopulated compared to that of 
SDSS ellipticals. We cannot formally rule out the 
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possibility that the most massive ULIRGs have 
been missed due to sample incompleteness, small 
volume densities, or possibly shorter (than the 
average) burst timescales for the most massive 
sources. However, one source with high velocity 
dispersion, Mrk 273, has been observed in our 
sample (see also Hinz & Rieke 2006). Even for 
Mrk 273, the total mass of the system (which 
also depends on its effective radius) is insufficient 
to classify it as a giant elliptical (of stellar mass 
> 5 x 10 11 Mq). Moreover, the fundamental plane 
analysis for all ULIRGs in our sample seems to 
suggest that ULIRGs do not account for the for- 
mation of a significant fraction of the giant ellip- 
tical galaxies in the local Universe. To quantify 
this statement, we apply a KS test to the ULIRG 
remnants and the (literature-selected) giant Es we 
used to construct the fundamental plane, acknowl- 
edging the different space-volumes they occupy. 
We find that the probability that ULIRG rem- 
nants and the giant ellipticals originate from the 
same population is extremely small, 10 -9 (see also 
Genzel et al. 2001); this probability indicates that 
despite the small number statistics, the confidence 
of this result is very high. For their morphologi- 
cally selected remnant sample, Rothberg & Joseph 
(2006a) find a higher, but still insignificant, prob- 
ability of 10%. 

To further assess how typicially ULIRGs form 
giant Es, we used the mass function of local Es to 
estimate what fraction / g i an t of the sources with 
mass > 10 10 Mq (mass threshold similar to that 
of ULIRGs) correspond to giant ellipticals (of stel- 
lar mass > 5 x 10 11 Mq). Given that the fraction 
may significantly vary according to the adopted 
mass function, we used the best fit to the K- 
band mass function from various authors in the 
literature. Giant ellipticals comprise 18.7% of lo- 
cal Es with mass > 10 10 M Q for the mass func- 
tion that Bell et al. (2003) derived for Es in 
the Sloan Digital Sky (SDSS) and the Two Mi- 
cron All Sky (2MASS) surveys assuming a diet 9 
Salpcter initial mass function (IMF). The frac- 
tion / giant is 20.1% for the 2MASS if -band lu- 
minosity function of local Es (Kochanek et al. 
2001), which we convert to a mass function us- 
ing m/L = 1.32 Mq/ Lq (Cole et al. 2001). This 

'In this initial mass function, the number of low-mass stars 
is reduced leading to a total stellar mass reduced by 30%. 



if-band m/L ratio is computed for both early- and 
late- type galaxies in the 2MASS and 2dF Galaxy 
Redshift Surveys, assuming a Salpeter IMF. An 
estimate of the m/L ratio for local ellipticals only 
can be derived from the work of Cappellari et al. 
(2005a) using the mean i-band m/L value and I-K 
color of the SAURON project sources. Combining 
the 2MASS if -band luminosity function with this 
early-type-galaxy m/L value (0.75 Mq/Lq) yields 
/ g iant = 10.9%. Although /giant is between 10% and 
20%, we do not observe any source of m > 5 x 10 11 
Mq in our sample of 29 ULIRG remnants. This 
fact indicates that local giant ellipticals cannot be 
accounted for, if assumed to originate only from 
local ultraluminous mergers. 

One possible scenario for the formation of local 
giant Es is that most of these objects have formed 
at higher redshifts (e.g. Swinbank et al. 2006). 
This scenario is based on the fact that the more 
massive early-type galaxies are, the faster and 
the earlier their stellar populations have formed 
(Thomas et al. 2005; van Dokkum 2005 and 
references therein). It is reinforced by the argu- 
ment that local giant Es are mainly located in 
clusters whereas ULIRGs lie in the field (Tacconi 
et al. 2002). Another likely scenario, which may 
hold simultaneously, is that other types of mergers 
(e.g. E-S or E-E ones) account for the formation 
of ^60% of elliptical galaxies (Khochfar & Burk- 
crt 2003; Kaviraz et al. 2006) at low redshift. 
E-S mergers enhance the star-formation activity 
already existing in one of the merging components 
and, for adequate amounts of gas, may appear as 
LIRGs. Gas-depleted E-E mergers (also known as 
"dry" mergers) are not easy to detect observation- 
ally (e.g, van Dokkum 2005; Bell et al. 2006) 
since they do not have other unambiguous sig- 
natures such as luminous/ultraluminous phase(s) 
beyond their perturbed morphologies; therefore, 
they can mainly be traced by visual identifica- 
tion. Dry merging is now believed to be a key 
element in the formation of local giant Es (Naab 
et al. 2006a) and is a likely explanation for the 
locations of some optically selected mergers in the 
fundamental plane (see § 5.2). 

Rothberg & Joseph (2006a) recently compiled 
a sample of local merger remnants selected by vi- 
sual inspection. Such a sample probably comprises 
merger remnants of all possible (S-S; E-S; E-E) 
categories. Rothberg & Joseph (2006a) measured 
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a from Ca II triplet (optical) spectroscopy for 38 
remnants, of which, 10 were LIRGs and 2 were 
ULIRGs. They found that the hypothesis that the 
(U)LIRGs in their sample and giant ellipticals can 
be drawn from the same parent population has a 
probability of 10%. Lake & Dressier (1986) also 
obtained optical velocity dispersion measurements 
for 13 merger remnants. Of those, 10 have In- 
frared Astronomical Satellite (IRAS) fluxes; none 
is a ULIRG and only one source, AM 0921-225, 
is a LIRG 10 . Rothberg & Joseph (2006a) com- 
puted the probability that the remnants of Lake 
& Dressier (1986) and giant Es originate from the 
same parent population and found it to be 5%. 

The conclusion from the comparison of the 
optically-selected remnants of Rothberg & Joseph 
(2006a) and Lake & Dressier (1986) with the IR- 
selected remnants in our study (and those of Gen- 
zel et al. 2001 and Tacconi et al. 2002) is what 
is expected from a theoretical point of view: while 
mergers account for the formation of galaxies of 
various morphological classes and mass contents 
(e.g., Schweizer & Seitzer 1992; Springel & Hern- 
quist 2005; Springel et al. 2005a), ultralumi- 
nous IR mergers seem to have a specific output 
(moderate-mass ellipticals) originating from their 
specific input (mostly sub-m* spirals). 

5.4. Wavelength Dependence Of The Stel- 
lar Velocity Dispersion Measurement 

Systematic differences between the measured 
and the actual values of the host dispersion may 
arise when extracting a from different wave- 
lengths, i.e. from the Ca II triplet in the op- 
tical (e.g. Rothberg & Joseph 2006a) and the 
CO bandheads in the NIR. The host kinematics 
extracted from the NIR bandheads are often rep- 
resentative of young stellar populations. These 
populations could still be linked to the gas from 
which they formed (which is believed to settle 
into a disk earlier than the stars in the progen- 
itor disks; Mihos & Hernquist 1996), and have 
less perturbed orbits than the old stars. Host dis- 
persions extracted from the CO bandheads could 
be systematically lower than those of the merg- 

'We computed the IR luminosity of these sources with the 
formula of Sanders & Mirabel (1996) using mid-infrared 
fluxes from Moshir et al. (1990), The Point Source Catalog 
(1988), and the NASA/IPAC Extragalactic Database 



ing bulges. However, the low V TO t/a ratio that 
we measure for our ULIRGs is a good indication 
that even the young stellar populations are signifi- 
cantly heated, and therefore, systematics originat- 
ing from NIR dispersion measurement are unlikely 
to have a major effect on our conclusions. 

In the case of the Ca triplet, systematic errors 
are mainly related to the presence of dust. Due 
to scattering of light from dust particles, photons 
originating from high-velocity stars in the center 
of the galaxy will be scattered into random lines of 
sight, biasing the LOS velocity distribution (Baes 
& Dejonghe 2002). The effects of dust in the 
observed stellar kinematics depend upon the dust 
mass and extent (Baes & Dejonghe 2002). Since 
ULIRGs are highly obscured systems, the stellar 
dispersions extracted from the Ca triplet may be 
systematically higher than their intrinsic values. 

Such possible systematics have not been inves- 
tigated in ULIRGs. Silge & Gebhardt (2003) have 
attempted to quantify the discrepancy between 
dispersion estimates derived from the NIR and the 
optical regimes for local quiescent galaxies. While 
they found no significant difference for elliptical 
galaxies, the systematics were non-negligible in SO 
lenticulars, with the largest difference in a (~30- 
40%) found in the most massive sources of their 
sample. The opposite conclusion is found from 
the velocity dispersion of ellipticals and Seyfert- 
type AGNs in Oliva et al. (1995) and Oliva et al. 
(1999). The a values for the objects in their sam- 
ple were systematically greater in the NIR than 
in the optical. The difference was large for type 

1 Seyferts (up to a factor of 2) but low for type 

2 Seyferts and ellipticals (<10%). However, one 
of the uncertainties often underestimated in these 
comparisons is that the spectral resolution in the 
NIR was significantly lower than that in the op- 
tical. For Silge & Gebhardt (2003), R was 2300 
in the if-band; for Oliva et al. (1995; 1999) it 
was equal to 1900 and 2500 in the H— and in the 
K— band respectively. To get a better estimate 
of the possible systematic errors in er, we com- 
pare the results for LIRGs with a Ca II triplet 
a measurement (Rothberg & Joseph 2006a) that 
also have CO observations of relatively high spec- 
tral resolution (e.g., equal or better than ours). 
These sources (NGC 1614, NGC 2623, NGC 4194, 
Arp 193, and IRAS 20551-4250) have been ob- 
served in the NIR by James et al. (1999), Gen- 
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zel et al. (2001), Rothberg & Joseph (2006a), 
and more recently, Hinz & Rieke (2006). For 
most sources, the optical dispersions were higher 
than the NIR ones, but the systematic difference 
was small (7% on average). On the other hand, 
significant deviations (of a factor of 2) exist be- 
tween the NIR results of James et al. (1999) and 
Hinz & Rieke (2006) and the older (also NIR but 
lower-resolution) data of Shier & Fischer (1998) 
for NGC 1614 and NGC 2623. It is therefore not 
clear whether the measurement of a from the dif- 
ferent features would systematically yield differ- 
ent results for ULIRGs, and if yes, whether they 
would reflect real dynamic trends instead of mea- 
surement errors. Errors may also be easily intro- 
duced from the use of different stellar templates 
by various authors. 

We investigate whether the use of the Ca II 
triplet for the measurement of a may affect our 
results on the formation of (giant) Es under the 
hypothesis that the optical a values can be in- 
deed greater than the NIR ones due to population 
and extinction effects. We compute the maximum 
possible increase that can be applied to the dis- 
persions of the sources in our sample, which is the 
maximum reported deviation (40%) in the results 
of Silge & Gebhardt (2003), corresponding to an 
addition of ~60 kms -1 in a or a shift of 0.15 in 
the (logarithmic) horizontal axis of the fundamen- 
tal plane. It is clear that this shift brings merger 
remnants closer to, but not into, the region popu- 
lated by giant ellipticals. We also recompute the 
masses of the individual remnants in our sample 
after increasing their dispersion by 40%. The av- 
erage mass is now 1.69xl0 n M , which is again 
close to m* and there is still no source more mas- 
sive than 5x10 11 M Q . Therefore, our main con- 
clusions are insensitive to the wavelength depen- 
dence of the stellar velocity dispersion measure- 
ment. Furthermore, the Kot(obs)/<7 ratios 11 of 
SO lenticulars are typically greater than those of 
ULIRGs; e.g, for the combined sample of Pizzcla et 
al. (2005) and Cappellari et al. (2005b) the ratio 
is twice that of the remnants in our sample (0.6) 
indicating that deviations of the order 30-40% be- 
tween the CO and Ca dispersion measurement are 
probably high for ULIRGs. 

not corrected for inclination effects 



6. Black Holes in ULIRGs 

6.1. A Picture Of The M bh -ct Relation 
Time Evolution 

The size of a black hole seems to be closely 
linked to the depth of the potential well in which 
it forms and grows (Ferrarese & Merritt 2000). 
This is reflected, e.g., in the local Mbh — & rela- 
tion, the correlation between the black hole mass 
and the stellar dispersion in the bulge of the host 
galaxy (Ferrarese & Merritt 2000; Gebhardt et 
al. 2000). Since the Mbh ~ o relation is found 
for virializcd bulges, it should be valid at the end 
of the merger process, when the bulge stellar kine- 
matics have reached their final dynamical state 
and AGN winds and supernovae feedback have 
expelled the gas away from the nucleus, prevent- 
ing further BH growth, terminating the starburst 
phase, and making the system resemble an ellipti- 
cal galaxy. The merger remnants in recent simula- 
tions (Di Matteo et al. 2005; Springel et al. 2005; 
Robertson et al. 2005) are able to reproduce the 
Mbh — cr relation by subjecting a portion of their 
interstellar gas to accretion and feedback. 

For ULIRGs however, the conversion of the host 
velocity dispersion into Mbh carries the uncer- 
tainty of applying the Mbh — cr relation to sys- 
tems out of dynamical equilibrium. It is not yet 
known if, or to what extent, the Mbh — o relation 
is valid between the first encounter and shortly 
(< 10 8 yrs) after the final coalescence, when most 
of the ultraluminous infrared activity occurs. 

To investigate whether (and under which condi- 
tions) merging disk galaxies fall on the A'/bh — c 
relation during interaction phases prior to coales- 
cence, we ran gas-rich merger simulations (details 
presented in Naab et al. 2006b) that have al- 
ready been discussed briefly in Paper I. To in- 
clude the effects of a dissipative component we re- 
placed 10% of the stellar mass in the initial disks 
with isothermal gas at a temperature of approx- 
imately 10 4 K. The initial scale length of the gas 
disk was equal to that of the stellar disk, h. Each 
galaxy had a stellar bulge with 1/3 of the disk 
mass and was embedded in a pseudo-isothermal 
halo to guarantee a flat rotation curve at large 
radii. All galaxies approached each other on nearly 
parabolic orbits with a pericenter distance of two 
disk scale lengths. The spin orientations were both 
prograde and retrograde. The evolution of the 
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stars and the gas was computed with the N-body 
smoothed-particle-hydrodynamics code VINE us- 
ing an isothermal equation of state for the gas. For 
this study we analyzed mergers with 16 different 
initial disk orientations (geometries 1-16 of Naab 
& Burkert 2003) and mass ratios 1:1 and 3:1. 

We analyzed snapshots in the orbital plane ap- 
proximately every half-mass rotation period of the 
more massive disk to follow the time evolution 
of the projected central stellar velocity dispersion 
and the gas accretion history onto the center of 
the system. At each snapshot, we calculated Mbh 
from the velocity dispersion assuming that the 
-Mbh — cr relation is always valid. As in Paper 
I, we used the Tremaine et al. (2002) formula 

M BH = 1-35 x 10 8 [cr/200] 4 - 02 M Q , (5) 

where a is in units of kms -1 , to calculate the 
black hole masses. To quantify the gas accretion, 
we computed the total gas mass M gas that has 
ever reached a radius of O.l/i (i.e. twice the reso- 
lution element of the simulations) relative to the 
mass center of the system. By calculating the ratio 
between the quantities Mbh and M gas we derive 
the accretion efficiency e that is needed to main- 
tain the black hole mass on the Tremaine et al. 
(2002) formula at any time. We define accretion 
efficiency e as the fraction of the gas that accretes 
onto the black hole from the gas accumulated into 
the nuclear region. 

In Fig. 9 we show the evolution of the ac- 
cretion efficiency (needed to keep the black hole 
on the Mbh — cr relation) as a function of time 
(left column) and nuclear separation (right col- 
umn) . Time is given in units of the half- mass rota- 
tion period, Tx/2, of the more massive progenitor 
disk. The nuclear separation is given in disk scale 
lengths of the more massive disk. The models are 
scale- free; scaling the parameters to e.g., those of 
the Milky Way corresponds to a unit time of 13 
Myrs and a unit length of 3.5 kpc. Scaling them 
to the mean scale length and rotational velocity 
of all individual progenitors in our binary ULIRG 
sample yields a unit time of 13 Myrs and a unit 
length of 1.3 kpc. The efficiency is averaged over 
the 16 initial geometries for 1:1 and 3:1 mergers 
(upper and lower panels). The diagonally shaded 
area indicates the spread due to the varying initial 
disk geometries. 

According to the left panels of Fig. 9, before 



nuclear coalescence, the estimated black hole mass 
and the gas accumulated at the center of the sim- 
ulation already scale linearly. In other words, if 
we assume that the Mbh — cr relation is valid at 
any time during the merger, then we find that the 
accretion efficiency remains constant shortly after 
the first encounter (shown in vertical dashes) until 
relaxation. Vice-versa, if we assume that e stays 
constant during the merger and require only that 
the Mbh — cr relation be valid at relaxation 12 , we 
find that the black hole mass can be computed 
from the Mbh — cr relation from midway of the 
first encounter through relaxation. 

The implication of these simulations is that the 
Mbh — cr relation can be used to compute black 
hole masses in merging systems, from mid- way be- 
tween first encounter and nuclear coalescence until 
the remnant becomes an elliptical galaxy, as long 
as e remains constant. This is true for both cqual- 
(1:1) and unequal- (3:1) mass major mergers. A 
constant efficiency e conceptually corresponds to 
subjecting a specific gas fraction to AGN feeding, 
an assumption also made for models that include 
star formation and ISM feedback (e.g. Springel 
et al. 2005). However, it is likely that on short 
timescales (e.g. during specific accretion events) 
the efficiency e may vary significantly. Therefore 
individual black hole estimates calculated from the 
Mbh — cr relation carry this uncertainty. 

In the right panels of Fig. 9, the "folding" of the 
efficiency e at large nuclear separations is due to 
the fact that the galaxies move towards apogalac- 
ticon before they fall back together. Clearly, at 
these merger phases, the use of the Mbh — cr rela- 
tion is misleading. For a constant e, the Mbh — cr 
relation begins to be valid beyond 5 disk scale 
lengths, which corresponds to a nuclear separation 
^7 kpc for our ULIRGs (see Paper I). If the as- 
sumption of a constant e holds, the black hole esti- 
mates for ~ 2/3 of the binary sources may be con- 
sidered reasonable (but still lower limits on their 
final values). For the ULIRG remnants, comput- 
ing BH masses from stellar dispersions seems also 
plausible. This is further supported by the fact 
that the kinematics of the latter are expected to 
have almost settled to their relaxation values (Mi- 
hos 2000). 

This assumption is reasonable since ULIRGs will eventually 
form Es. 
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6.2. Black Hole Sizes and Accretion Rates 

We calculate the individual black hole masses of 
the sources in our sample using the Tremaine et al. 
(2002) formula and tabulate the results in Table 3. 
The mean black hole mass of the ULIRG remnants 
is 8.4 x 10 7 Mq. Since our sources are still in an ul- 
traluminous IR emission phase, gas and dust are 
still present in the nuclear region and will con- 
tinue accreting onto the black hole. However, be- 
yond coalescence, the timescales for further gas 
accretion onto the black hole are probably shorter 
than those of the pre-coalescence starbursts (e.g. 
Springel et al. 2005). Therefore, the black hole 
masses of the fully relaxed systems will be some- 
what larger than those presented in Table 3. Sim- 
ple gas content (typical gas mass of ULIRGs is 
5 x 10 9 M & ; Downes & Solomon 1998) and ac- 
cretion efficiency (0.01, see § 6.1) arguments indi- 
cate that the additional increase of the remnant 
BH masses will not exceed 10 s M Q by the time 
the ultraluminous activity ends. In this order-of- 
magnitude calculation we have ignored the fact 
that only part of the gas reaches (or remains in) 
the center of the merging system, further reducing 
the upper limit on the final black hole mass. 

We calculate the BH mass each source would 
have, if it were accreting at the Eddington rate 

L E dd/L Q = 3.8 x 10 4 M Ed d/M Q . (6) 

To calculate the Eddington BH mass M^dd, we 
assign to the Eddington luminosity ^Edd half of 
the luminosity emitted in the IR, given that some 
ULIRGs are largely starburst- while others are 
AGN- powered (Genzel et al. 1998; Rigopoulou 
et al. 1999; Veilleux et al. 1999; Joseph 1999; 
Sanders 1999). The Eddington efficiency 7yEdd, 
the ratio of the Eddington Mbh estimate over the 
dynamical black hole mass, is given in Table 3. 
The mean Eddington efficiency of the merged- 
ULIRG sample is 0.89, with a wide range of val- 
ues. While statistically appropriate, the assump- 
tion that 50% of the IR luminosity originates from 
the AGN may make some sources appear as if ac- 
creting at super-Eddington rates (see Table 3). If 
we exclude these sources from the statistics and 
recalculate the average efficiency, we find it to be 
0.39. 

The inferred accretion rates may be higher than 
in reality if an overestimated fraction of the IR 



emission is assigned to the AGN. This could occur 
if, for example, the starburst is still the dominant 
source of the IR luminosity after the merging nu- 
clei coalesce. The relative strengths of the star- 
burst (Lsb) and the AGN (Lagn) luminosity at 
each merger phase are uncertain. However, the 
MIR ISO spectroscopic study of local ULIRGs by 
Genzel et al. (1998) and Rigopoulou et al. (1999) 
indicates that most ULIRGs are starburst domi- 
nated systems, implying that probably less than 
50% of the IR luminosity should be assigned to the 
AGN. Veilleux et al. (2002) have found that the 
strength of nuclear continuum emission increases 
with decreasing nuclear separation. Therefore, the 
Eddington efficiencies of the progenitors presented 
in Paper I may be systematically overestimated 
compared to those of the remnants. Future papers 
presenting SPITZER MIR spectroscopy of local 
ULIRGs will indicate the appropriate luminosity 
fraction that needs to be assigned to the AGN for 
ULIRGs before and after nuclear coalescence. 

7. Conclusions 

We have acquired NIR spectroscopic, long-slit, 
data of 54 ULIRGs at a variety of merger phases to 
trace the evolution of their host dynamical prop- 
erties. From the analysis of the stellar kinematics 
in 29 ULIRG remnants, we find that: 

1. Indications of an increase of the stellar 
random motions exist as the merger ad- 
vances. The mean stellar <r, as mea- 
sured from the CO rovibrational bandheads, 
equals 142 kins" 1 for the binary sources 
(of mean nuclear separation 9.4 kpc, includ- 
ing IRAS 00456-2904 and IRAS 09111-1007) 
and 161 kms -1 for the remnants (of nuclear 
separation <1.5 kpc). This difference in 
the means of the pre- and post- coalescence 
distributions is marginally significant and 
requires more data to be accurately con- 
strained. This increase of the stellar dis- 
persion observed in ULIRGs corresponds to 
only a part of the dynamical heating that 
occurs during the merger, since the merger 
timescales are longer than those of the ul- 
traluminous starburst. 

2. The dynamical and structural properties of 
the remnants indicate that they originate 
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from mergers mainly of 1:1 and 2:1 progen- 
itor mass ratios. This confirms what we 
found in Paper I by directly measuring the 
masses of the individual progenitors of bi- 
nary ULIRGs. 

3. Ultraluminous mergers will mainly lead to 
the formation of moderate mass ellipticals 
(of stellar mass ~10 n M©). Depending on 
the definition of to*, ULIRGs are between 
sub- and ~ to* . ULIRGs are located in a re- 
gion of the fundamental plane of early-type 
galaxies different from that of local giant Es, 
indicating a different formation history for 
most of the latter. Local mass and disper- 
sion functions support this argument; how- 
ever, the current statistics cannot exclude 
the formation of few giant Es from ULIRGs. 

4. We have performed simulations to investi- 
gate whether a black hole mass-host disper- 
sion relation can be used to calculate the 
black hole masses of our ULIRGs. We find 
that already before nuclear coalescence, the 
mass of the gas that falls into the center of 
the merging system scales linearly with the 
black hole mass predicted by the Mbh — o" 
relation. However, this is only true (and con- 
versely) if the efficiency of gas accretion onto 
the BH from its surroundings remains con- 
stant with time. 

5. The black hole masses of the merged ULIRGs 
are of order 10 7 -10 8 M© and their accretion 
rates are high (Eddington efficiencies often 
> 0.5). If the AGN luminosity output of 
a ULIRG nucleus increases with time, our 
accretion rates in pre-coalescence ULIRGs 
may be overestimated relative to those in 
post-coalescence ULIRGs. 
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A. APPENDIX: Notes on individual sources 

IRAS 00091-0738: This source has two overlapping nuclei separated by 1.4 kpc. Archival Wide Field 
Planetary Camera 2 data of this source (PI Borne) also indicate an optical nuclear separation of 1.5 kpc. 

IRAS 00456-2904: According to the imaging analysis of Kim et al. (2002) and Vcilleux et al. (2006) this 
source is probably a binary system at projected nuclear separation of 22.2 kpc (converted to our cosmology). 
The redshift of the NE nucleus is not spectroscopically confirmed but its ambiguous morphology indicates 
an interaction. The data presented in this paper are for the more luminous (late-type-host) SW nucleus. 

IRAS 09111-1007: Due et al. (1997) find that this source has two widely separated interacting nuclei; the 
ultraluminous IR component is the western source. The nuclear separation measured from our acquisition 
images equals 40.4 kpc. 

IRAS 14348-1447: According to the NICMOS imaging of Scoville et al. (2000), this source has a projected 
nuclear separation of 6.0 kpc (converted to our cosmology). 

IRAS 14378-3651: Two components separated by 1.5 kpc are visible in the acquisition image of this 
source. 

IRAS 15250+3609: The two components of this source are separated by 0.8 kpc (Scoville ct al. 2000). 
IRAS 17208-0014: Hinz & Rickc (2006) measure the velocity dispersion of this source to be 125 ± 28 
kms' 1 . 

IRAS 20551-4250: Rothberg & Joseph (2006a) find a dispersion of 185 ±6 kms" 1 for this source. 

IRAS 23578-5307: Our acquisition image indicates the presence of two nuclei separated by 1.4 kpc. 

Arp 220: This source has a projected nuclear separation of 0.3 kpc (Scoville et al. 1998). 

NGC 6240: The two nuclei of this source are separated by 1.4 kpc and still show some discrete velocity 
patterns. The velocity dispersion of the system peaks between them; at the position of this peak the stellar 
kinematics may reflect localized motions of self-gravitating gas (Tecza et al. 2000; Genzel et al. 2001). For 
the velocity dispersion of this source we use the average of the values at two nuclei and at the internuclear 
peak (229 kms -1 ), which is close to the luminosity-averaged value, 225 kms -1 (Tecza et al. 2000). 

Mrk 273: This system has two nuclei separated by 0.8 kpc (Scoville ct al. 2000). Hinz & Rieke (2006) 
find a stellar velocity dispersion of 232 ±43 kms -1 for the north component of this source. 
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0.037 


12.13 


15,95 


40,40 


R 


13451 + 1232 b 


13:47:33 


12:17:23 


0.122 


12.28 


104,13 


80,120 


B 


14070 + 0525 


14:09:31.3 


+05:11:31 


0.264 


12.76 


-1,89 


120,120 


R 


14348-1447 c 


14:37:38.3 


-15:00:23 


0.0823 


12.3 


30 


240 


B 


14378-3651 c 


14:40:58.9 


-37:04:33 


0.068 


12.24 


-45 


80 


R 


15130-1958 


15:15:55.2 


-20:09:17 


0.109 


12.09 


-171,-81 


80,110 


R 


15250 + 3609 c 


15:26:59.4 


-35:58:38 


0.055 


11.99 


45, -45 


40,40 


R 


15327+2340 (Arp 220) c 


15:34:57.1 


23:30:11 


0.0181 


12.1 


52, -91 


120,120 


R 


15462-0450 


15:48:56.8 


-04:59:34 


0.100 


12.16 


179,-91 


180,160 


R 


16156+0146 b 


16:18:08 


01:39:21 


0.132 


12.04 


-50,-51,40,40 


60,60,60,60 


B 


16300 + 1558 b 


16:32:20 


15:51:49 


0.242 


12.63 


-1,89 


150,90 


B 


16504+0228 (NGC 6240) c 


16:52:58.9 


02:24:03 


0.0245 


11.8 


-158,-31 


20,20 


R 


17208-0014 c 


17:23:21.9 


-00:17:00 


0.0428 


12.33 


90,120 


30, 30 


R 


19254-7245 b 


19:31:21.4 


-72:39:18 


0.0617 


12.00 


-13,77 


60,60 


B 


20046-0623 b 


20:07:19.3 


-06:14:26 


0.0844 


11.97 


69,159 


60,60 


B 


20087-0308 c 


20:11:23.2 


-02:59:54 


0.106 


12.40 


-45,45 


40,40 


R 


20414-1651 


20:44:18.2 


-16:40:16 


0.087 


12.26 


54,144 


100,80 


R 


20551-4250 c 


20:58:26.9 


-42:39:06 


0.0428 


11.98 


-45,45 


60,60 


R 


21130-4446 b 


21:16:18.5 


-44:33:38 


0.0926 


12.02 


33 


40 


B 


21208-0519 b 


21:23:29 


-05:06:59 


0.13 


12.01 


-164,109,109 


60,60,60 


B 


21219-1757 


21:24:41.6 


-17:44:46 


0.112 


12.06 


-1,89 


50,40 


R 


21329-2346 b 


21:35:45 


-23:32:36 


0.125 


12.09 


31 


60 


B 


21504-0628 


21:53:05.5 


-06:14:50 


0.078 


11.92 


-39,59 


60,60 


R 


22491-1808 b 


22:51:49.2 


-17:52:23 


0.0778 


12.09 


-76,13,13 


60,60,60 


B 


23128-5919 b 


23:15:46.8 


-59:03:15 


0.045 


11.96 


-5,84,84 


40,40,40 


B 


23230-6926 


23:26:03.6 


-69:10:19 


0.106 


12.17 


-1,89 


60,60 


R 


23234+0946 b 


23:25:56.2 


10:02:50 


0.128 


12.05 


-64,25 


60,60 


B 


23365 + 3604 c 


23:39:01.3 


+36:21:10 


0.0645 


12.09 


45, -30 


15, 40 


R 


23578-5307 c 


00:00:23.6 


-52:50:28 


0.125 


12.10 


107,14 


60, 40 


R 



NOTE. — The coordinates, the redshift, the infrared luminosity, as well as the position angles and respective integration time for all our 
sources are presented in this Table. The data arc presented in this paper unless otherwise noted. 

a Pre- or post- coalescence classification, according to the scheme described in § 2.2. "B" refers to a binary (pre-coalescence) source. 
"R" refers to a remnant (in post-coalescence). Remnants may also have two nuclei, separated by no more than 1.5 kpc. "U" indicates 
that the classification is uncertain. 

b Source presented in Paper I. 

c Source presented in Genzcl et al. (2001) and Tacconi ct al. (2002). 
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Table 2 

Derived structural and photometric parameters 



Galaxy 




cllipticity 


A* 
<Pa 


-Tt c ff^-f\ 0(111(1) 


Mgai (^ bancT) 






(kpcj 




\ ) 


[kpc) 


mag ( ) 


mag ( ) 


00091-0738 


2.47(±0.21) 


0.220 


15 






16.8 h 


00262+4251 a 








3.4 (±1.0) c 






00397-1312 b 


2.04 (±0.76) 


0.35 


-25 




16.4 


15.2 h 


F01004-2237 b 


0.40 (±0.07) 


0.02 


26 




14.3 


12.9 k 


01388-4618 


1.62 (±0.03) 


0.074 


55 






14.7 I ' h 


01572+0009 


1.31 (±0.10) 


0.140 


81 


3.16 (±1.42) c 


14.6 d 


13.7 k 


F02021-2103 b 


5.38 (±3.95) 


0.34 


1 




17.1 


17.2 h 


04103-2838 b 


1.61 (±0.12) 


0.19 


-80 




15.6 


15.0 h 


04313-1649 b 


4.04 (±0.89) 


0.26 


86 


5.7 f 


18.5 


18.0 h 


05189-2524 b 


0.57 (±0.08) 


0.06 


77 


0.79 (±0.02) c 


13.6 


13.5 k 


09039+0503 b 


1.62 (±0.90) 


0.15 


-18 




16.2 


15.5 h 


11223-1244 


3.83 (±0.38) 


0.101 


70 






16.1 h 


12540+5708 b 


1.40 (±0.21) 


0.08 


-11 


0.2 ° 


13.4 


13.2 k 


13428+5608 c 


1.03 (±0.29) 


0.498 


-77 


1.15(±0.24) c 


14.5 d > h 


14.1 d ' h 


14070+0525 b 


3.62 (±0.53) 


0.17 


-19 




17.2 


16.2 h 


14378—3651 


0.36 (±0.02) 


0.031 


20 


0.67 e 


14.6 ' 


13.4 s ' h 


15130—1958 b 


1.62 (±0.27) 


0.23 


-88 




16.1 


15.0 h 


15250+3609 c 


2.10 (±0.09) 


0.207 


-72 


1.49(±0.12) c 


16.9 d ' h 


16.5 d ' h 


15462-0450 b 


5.57 (±1.38) 


0.048 


62 


1.6 f 


17.5 


16.2 h 


17208-0014 


1.69 (±0.08) 


0.196 


54 


1.63(±0.07) c 


15.9 d ' h 


15.6 d > h 


20087-0308 


1.87 (±0.22) 


0.324 


-86 






15.2 j' h 


20414-1651 b 


1.37 (±0.47) 


0.67 


3 




16.2 


15.7 h 


20551-4250 


1.32 (±0.25) 


0.113 


84 


3.3s 




16.9 j' h 


21219-1757 b 


4.19 (±3.29) 


0.14 


-50 




16.9 


15.6 h 


21504-0628 


1.95 (±0.26) 


0.164 


12 








23230-6926 


2.03 (±0.25) 


0.177 


41 




16.9 1 


16.2 j' h 


23365+3604 a 


4.8 (±1.0) c 












23578-5307 


3.96 (±1.36) 


0.447 


-80 








Arp 220 d 


3.00 






0.69(±0.03) c 


17.4 d - h 


16.9 d ' h 


NGC 6240 d 


1.44 






0.92(±0.23) c 


14.3 d ' h 


13.8 d > h 


00456-2904(SW) 


2.09 (±0.15) 


0.067 


-88 




16.3 


14.9 


09111-1007(W) 


2.18 (±0.36) 


0.397 


-26 








14348-1447(NE) 


3.03 (±0.50) 


0.237 


4 


3(±l) c 


18.3 j 


17.5 j 


14348-1447(SW) 


2.78 (±0.03) 


0.113 


7 


3(±l) c 


17.8 J 


16.8 j 



Note. — Structural and photometric properties of (mainly) ULIRG remnants. In this table we present 
data for newly observed sources during our VLT large program and we summarize observations carried out 
during our older normal programs (Gcnzel ct al. 2001; Tacconi et al. 2002) that are used in this work. The 
sources in this Table together with those presented in Paper I summarize the 54 sources in our sample. The 
H-band structural parameters R c g, ellipticity, and <j> a of the sources are derived from the ISAAC acquisition 
images unless otherwise noted in the first column. The angle <f> a is tabulated here only for the first slit; for 
the second slit it can be computed using the angle between the two slit positions given in Table 1. The galaxy 
surface brightness within the effective radius, /i ga i is derived from Veilleux ct al. (2006) and Kim ct al. 
(2002) respectively, unless otherwise noted. The H- and ii"-band /i ga i data arc presented after subtraction of 
the nuclear PSF for both bands; we indicate when individual PSF and galaxy magnitudes are not available in 
the literature. All data presented in this Table are derived from maximum-aperture photometry (enclosing 
the whole galaxy) and are not corrected for extinction effects. 

a For the sources observed using NIRSPEC and the slit monitoring camera, the structural parameters are 
not extracted; the slit projection drawn on the acquisition image does not allow for photometric analysis. 
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b Data taken from H-band NICMOS imaging by Veilleux et al. (2006). All quantities arc PSF-subtracted 
and are converted to our cosmology. 

Structural parameters extracted from NICMOS imaging data, kindly provided to us by Nick Scoville. 

d Based on Scoville ct al. (2000). 

°Effcctivc radius from Genzel et al. (2001) or Tacconi et al. (2002). 
f Data from the 1 Jy sample analysis (Veilleux et al. 2002). 
g Effective radius from Rothbcrg & Joseph (2004). 

h For these sources the relative strength of the PSF and the Sersic component fit to the underlying 
galaxy are not known. We assume that L ga i(H band) = 0.64Ltot{H band) and that L ga i(K band) = 
0A8L to t(K band), (Colina et al. 2001; Suracc & Sanders 1999; Veilleux et al. 2006). 

'/f-band magnitude taken from Colina et al. (2001). 

^magnitude taken from Due ct al. (1997). 

k i\-band magnitude taken from Surace & Sanders (1999). 

'ii"-band magnitude taken from Rigopoulou ct al. (1999). 
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Table 3 

Stellar velocities and black hole masses 



Source a apcr. a apcr. a V rot (obs) b apor. c V rot (obs) / 1 M BH d M BH (Edd.) ° r, Edd 



(IRAS) 


(km.- 1 ) 


( 


') 


(kpc) 


(km a" 1 ) 


(") 






(M ) 




( Af ) 




00091-0738 


131 (± 39) 





45 


1.1 








2 


46 


X 


10 7 


2 


04 


X 


10 7 


0.83 


00262+4251 


170 (± 15) 








< 15 




< 0.09 


7 


02 


X 


10 7 


1 


32 


X 


10 7 


0.19 


00397-1312 


106 (± 26) 





36 


2.0 


49 (± 17) 


0.58 


0.46 


1 


05 




10 7 


1 


05 




10 s 


9.94 


F01004-2237 


132 (± 29) 





29 


0.7 


22 (± 13) 


0.60 


0.17 


2 


54 




10 7 


2 


29 




10 7 


0.90 


01388-4618 


144 (± 10) 








130 (± 15) 




0.90 


3 


60 




10 7 


1 


32 


X 


10 7 


0.37 


01572 + 0009 


200 (± 60) 














1 


35 


X 


10 8 


4 


46 


X 


10 7 


0.33 


F02021-2103 


143 (± 21) 





37 


0.9 


42 (± 10) 


0.76 


0.29 


3 


50 


X 


10 7 


1 


35 


X 


10 7 


0.38 


04103-2838 


129 (± 40) 





29 


0.7 


4 (± 5) 


0.31 


0.03 


2 


32 


X 


10 7 


1 


86 


X 


10 7 


0.80 


04313-1649 


157 (± 21) 





38 


2.1 


31 (± 27) 


0.57 


0.20 


5 


10 


X 


10 7 


4 


67 


X 


10 7 


0.92 


05189-2524 


137 (± 16) 





29 


0.3 


70 (± 14) 


0.51 


0.51 


2 


95 


X 


10 7 


1 


62 


X 


10 7 


0.55 


09039+0503 


183 (± 38) 





58 


1.5 








9 


45 


X 


10 7 


1 


55 




10 7 


0.16 


11223-1244 


149 (± 27) 





73 


3.0 








4 


13 


X 


10 7 


5 


12 


X 


10 7 


1.24 


12540+5708 


120 (± 10) 








25 (±10) 




0.21 


1 


73 


X 


10 7 


4 


16 


X 


10 7 


2.40 


13428+5608 


285 (± 30) 








110 (±20) 




0.39 


5 


61 


X 


10 8 


1 


66 


X 


10 7 


0.03 


14070+0525 


139 (± 21) 





70 


3.8 


54 (± 19) 


0.77 


0.39 


3 


13 


X 


10 7 


7 


57 


X 


10 7 


2.42 


14378-3651 


153 (± 10) 








15 (±10) 




0.10 


4 


60 




10 7 


1 


66 




10 7 


0.36 


15130-1958 


177 (± 39) 





58 


1.3 


33 (± 21) 


0.66 


0.19 


8 


26 


X 


10 7 


1 


62 


X 


10 7 


0.20 


15250+3609 


150 (± 10) 








60 (±15) 




0.40 


4 


25 


X 


10 7 


1 


32 


X 


10 7 


0.31 


15462-0450 


169 (± 38) 





58 


1.2 








6 


86 


X 


10 7 


1 


90 


X 


10 7 


0.28 


17208-0014 


229 (± 15) 








110 (± 20) 




0.48 


2 


33 


X 


10 s 


2 


63 


X 


10 7 


0.11 


20087-0308 


219 (± 14) 








50 (± 15) 




0.23 


1 


94 




10 8 


3 


31 




10 7 


0.17 


20414-1651 


187 (± 32) 





44 


0.8 


96 (± 38) 


0.88 


0.51 


1 


03 




10 8 


1 


82 




10 7 


0.18 


20551-4250 


140 (± 15) 








40 (± 10) 




0.29 


3 


22 


X 


10 7 


1 


32 


X 


10 7 


0.41 


21219-1757 


121 (± 11) 





58 


1.3 








1 


79 


X 


10 7 


1 


51 


X 


10 7 


0.84 


21504-0628 


90 (± 31) 





32 


0.5 


9 (± 28) 


0.48 


0.10 


5 


45 


X 


10 6 


8 


30 


X 


10 6 


1.52 


23230-6926 


143 (± 14) 





48 


1.1 


23 (± 13) 


0.73 


0.16 


3 


50 


X 


10 7 


1 


44 


X 


10 7 


0.41 


23365+3604 


145 (± 15) 








< 15 




< 0.10 


3 


71 


X 


10 7 


1 


62 


X 


10 7 


0.44 


23578-5307 


190 (± 70) 














1 


10 




10 8 


1 


66 




10 7 


0.15 


Arp 220 


164 (± 10) 








185 (± 30) 8 




0.84 


6 


08 




10 7 


1 


66 




10 7 


0.27 


NGC 6240 


229 (± 43) h 








240 (± 108) h 




1.05 


2 


33 


X 


10 8 


8 


30 


X 


10 6 


0.04 


00456-2904(SW) 


162 (± 25) 








45 (± 10) 




0.28 


5 


79 


X 


10 7 


1 


73 


X 


10 7 


0.30 


09111-1007(W) 


112 (± 18) 





29 


0.3 


68 (± 16) 


1.17 


0.61 


1 


31 


X 


10 7 


9 


53 


X 


10 6 


0.73 


14348-1447(NE) 


170 (± 14) 








60 (± 20) 




0.35 


7 


02 


X 


10 7 


1 


92 


X 


10 7 


0.27 


14348-1447(SW) 


150 (± 25) 








50 (± 15) 




0.33 


4 


25 


X 


10 7 


3 


33 


X 


10 7 


0.78 



NOTE. — The stellar central velocity dispersion, rotational velocity, and the V ro ^/cr ratio of the ncwly-obscrved sources are derived from the 
spectra of Fig. 1 with the aid of the parameters of Table 2. The dynamical and Eddington black hole mass of each ULIRG and the ratio of 
the two arc also presented here. 

Radial extent of the central aperture used for the extraction of a. tabulated both in angular and linear distances. 

"The observed rotational velocity value presented in this column is corrected for angular deviations from the major axis of rotation but 
not for inclination effects. 

c The center of the outer aperture used for the extraction of V ro t . The annular width of this aperture is equal to the radius of the 
corresponding central aperture given in the third column. 

d Dynamical black hole masses estimated from their relation to the bulge dispersion (Trcmainc ct al. 2002). 

e Eddington black hole mass, calculated by attributing 50% of Lj R to the AGN. 

f Ratio of Eddington black hole mass over dynamical black hole mass. 

SWe use as rotational velocity of Arp 220 that of the east component, since that of the west component is only a lower limit (Genzel et 
al. 2001). 

h Data are taken from Tccza ct al. (2000). The average of the velocity dispersion values at the two nuclei and the internuclear region is 
tabulated here. This value is also close to the luminosity- weighted average of the two nuclei which equals 225 kms — 
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Fig. 1. — The iJ-band spectra of the ULIRG remnants from this study. The stellar template, convolved 
with Gaussians that represent the LOS broadening function of the sources, is overplotted as a solid line. All 
spectra are shifted to rest frame. 
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Fig. 1 continued. 
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Fig. 2. — The ISAAC -ff-band raw (left panel) and PSF-deconvolved (right panel) acquisition image of IRAS 
23578-5307. 
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Fig. 3. — Left panel: The distributions of stellar dispersions in ULIRGs before and after nuclear coalescence. 
Right panel: The ULIRG parent-population distribution used in our Monte Carlo simulations. 
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Fig. 4. — Ratio of observed rotational over dispersion velocity versus ellipticity for merger remnants (from 
Naab & Burkert 2003 and Naab et al. 2006b). Each panel corresponds to mergers of different progenitor 
mass ratios (cases from 1:1 to 4:1 are studied). The bold line, solid line, and dotted contours correspond 
respectively to the 90%, 70%, and 50% probability of finding a collisionless merger remnant in the enclosed 
region. The dark- and light-gray shaded areas indicate the region in the diagram where gas-rich merger 
remnants are expected to be found at 90% and 50% probability levels. The dashed line shows the theoretical 
values for an oblate isotropic rotator. ULIRG remnants are shown as triangles. 
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Fig. 5. — The R c ff-<r projection of the early-type galaxy fundamental plane. The data for the giant boxy and 
moderate-mass disky Es (squares and circles respectively) are taken from Bender et al. (1992) and Faber et 
al. (1997). More (cluster) Es (open circles) are from Pahre (1999). For viewing clarity, the various types of 
mergers are plotted in separate panels. The ULIRG remnants (29 from this study and 2 from Rothberg & 
Joseph 2006a) are plotted as triangles (left panel). LIRGs (diamonds) from Shier & Fischer (1998), James 
et al. (1999), Rothberg & Joseph (2006a), and Hinz & Rieke (2006) and other (visually-selected) mergers 
from Rothberg & Joseph (2006a) are plotted in diamonds and open-crossed diamonds respectively (right 
panel). The effective radii of all merger remnants used in this figure are equal to the averages of their NIR 
measurements, if more than one is available. 
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Fig. 6. — Upper panel: The 3-dimcnsional view of the if-band fundamental plane of early-type galaxies. 
Local early-type galaxies (mainly from the Coma and Virgo clusters) are plotted as open circles (Zibetti et 
al. 2002). 21 ULIRGs from this study are plotted as triangles. All data are PSF removed (see Table 2). No 
extinction corrections have been applied to the fluxes. The effective radii used are as in Fig. 5 Lower panel: 
The fundamental plane, viewed as in Pahre (1999). 
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Fig. 7. — Upper panel: The 3-dimcnsional view of the if-band fundamental plane of early-type galaxies. 
The symbols and effective radii used arc identical to those in Fig. 5. A PSF estimate has been removed to all 
photometric data initially taken from Kim et al. (2002) (see Table 2). No extinction corrections have been 
applied to the data. From the ULIRGs plotted in this figure, 25 are from this study and 2 from Rothberg & 
Joseph (2006a). Lower panel: The fundamental plane, viewed as in Pahre (1999). 
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Fig. 8. — The number density of sources as a function of their stellar velocity dispersion is plotted in this 
figure. The solid line corresponds to the number density per a of SDSS ellipticals, computed from their 
dispersion function (Sheth et al. 2003). The number density per a of ULIRG remnants, plotted as a 
histogram, is calculated by multiplying the % fraction of our remnants that resides in each a bin with the 
local volume density of ULIRGs from Sanders et al. (2003). To facilitate the comparison, the ULIRG 
histogram is further normalized so that its mean has the same number density as that of the elliptical n(a) 
distribution. 
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Fig. 9. — Evolution, during a disk-galaxy merger, of the average black-hole mass inferred from the central 
(line-of-sight) stellar velocity dispersion of the simulated galaxies. The black hole mass is in units of the total 
gas mass accreted onto the center of the simulations. In the left column, the evolution is given as a function 
of time. Time is plotted in units of the half-mass rotation period, 7i/ 2 , of the more massive progenitor 
disk. Open squares represent disks that are still separated, filled dots indicate fully merged systems. The 
spread originating from the initial disk orientations is indicated by the diagonally shaded area. The vertically 
shaded area indicates the first data point after the first encounter. In the right column, the evolution is 
plotted as a function of nuclear separation of the interacting galaxies. The nuclear separation unit is the disk 
scale length h of the more massive disk. We begin plotting data points for 1 step prior to first encounter, as 
in the left panels. We show 1:1 mergers in the upper and 3:1 mergers in the lower panels. 
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